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Recent lattice calculations of hadron structure functions are described. 



1. INTRODUCTION 

Since the late 1960's, deep-inclastic scatter- 
ing of unpolarised leptons and nucleons has been 
studied experimentally. The inclusive cross sec- 
tion for electroproduction can be described in 
terms of the structure functions Fi and F% re- 
lated to the densities of quarks and gluons in the 
nucleon (see, e.g., Ref. jlj). There is also some 
more indirect information on pion structure func- 
tions. Deep-inelastic scattering of longitudinally 
polarised leptons and nucleons allows one to mea- 
sure two more structure functions, g\ and g 2 - 

The observed scaling violations could be un- 
derstood in the framework of perturbative QCD 
The computation of the structure functions them- 
selves requires a nonperturbative method, e.g. lat- 
tice gauge theory and Monte Carlo simulations 
|2|. In this talk I shall describe recent attempts 
to calculate hadron structure functions on the lat- 
tice. An alternative method using the Hamilto- 
nian version of lattice gauge theory is discussed 
in Ref. §. 

Instead of calculating the structure functions 
directly, we compute hadronic matrix elements 
of certain operators which, through the operator 
product expansion, are related to moments of the 
structure functions. In the deep-inelastic limit, 
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the leading contribution is given by 
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Here n = 2,4, 6, . . .; / labels the contributing op- 
erators, and fi is the renormalisation scale needed 

to define the (reduced) matrix elements i>„ . The 
f f ) 

Wilson coefficients c\ ' are calculated in pertur- 
bation theory and depend on /i 2 /Q 2 as well as on 
the coupling constant g(fi). 

The operators we have to deal with are the "2- 
quark operators" (for f = q = u, d) 
n-i 
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(D = covariant derivative) and (for / = g) purely 
gluonic operators, which contribute only in the 
flavour-singlet sector. For the nucleon the matrix 
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elements Vn are then defined by 
traces \p, s) 
= 2v { n f) ( Plll '--p^- traces) , 
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where \p, s) denotes a nucleon state with momen- 
tum p and spin vector s (s 2 = — m 2 ). Sym- 
metrisation of all indices (indicated by {• • •}) and 
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subtraction of traces are needed to obtain oper- 
ators transforming irreducibly under the Lorentz 
group, i.e. operators of definite twist = 2. Anal- 
ogous results exist for pion and rho as well. 

( f) 

Within the parton model the u„ ' are inter- 
preted as average values of powers of the fraction 
of the hadron momentum carried by the parton 
(quark of flavour q or gluon): 
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For the polarised nucleon structure functions 
one obtains in the chiral limit: 
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for n = 0, 2, 4, . . . and 
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for n = 2, 4, . . . with Wilson coefficients e/^ and 
reduced matrix elements defined by [Q 
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— traces) . 
Here we need the operators (for / = q) 

= (5) 97.75 K w • • • K p „g . 

Note that dfl corresponds to twist 3. Within the 
parton model, |&q 9 ' = Aq is interpreted as the 
fraction of the nucleon spin carried by the quarks 
of flavour q. 

Now the challenge for lattice QCD is the cal- 
culation of the reduced matrix elements v n , a„, 
d n . So we have to compute forward hadron ma- 
trix elements of (relatively complicated) compos- 
ite operators. 



2. LATTICE CALCULATION 

These matrix elements are calculated from 
three-point correlation functions. With suitable 
operators B, B for the particle to be studied, 
e.g. the nucleon, we can write schematically for 
t > t > 



(B(t)O(r)B(0)) 

= (0\Be- H ^- T ^Oe- HT B\0) 

= (0\B\N)e- ENt (N\B\0)(N\O\N) 
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on a lattice with time extent T — > 00. Corre- 
spondingly we have for the two-point function 



(B(t)B(0)) = (0\Be- m B\0) 
' =°° {0\B\N)e~ ENt {N\B\0) 
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So we can determine the desired matrix elements 
from ratios of the form 
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(B(t)O(r)B(0)) 
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which should be independent of r for <C r <C t. 

The bare lattice operators have of course to 
be renormalised and may mix with other oper- 
ators in the process of renormalisation. In the 
cuclidcan continuum we should study operators 
like 
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or rather 0(4) irreducible multiplets with definite 
C-parity. In particular, we obtain twist-2 opera- 
tors by symmetrising the indices and subtracting 
the traces. In the flavour-nonsinglet case they do 
not mix and arc hence multiplicatively renormal- 
isable. 

Working with Wilson fermions (as we do) it 
is straightforward to write down lattice versions 
of the above operators. One simply replaces the 
continuum covariant derivative by its lattice ana- 
logue. However, 0(4) being restricted to its finite 
subgroup H(4) (the hypercubic group) on the lat- 
tice, the constraints imposed by space-time sym- 
metry are less stringent than in the continuum. 
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In particular, a multiplet of operators which is ir- 
reducible with respect to 0(4) will in general de- 
compose into several irreducible H(4) multiplets 
and the possibilities for mixing increase (^||. 

3. THE SIMULATIONS 

Our simulations were performed on Quadrics 
parallel computers. We have worked on lattices 
of size L 3 x T = 16 3 x 32 and 24 3 x 32 at = 6.0 in 
the quenched approximation. For more technical 
details see Refs. Most of our data were ob- 

tained with the Wilson action for the gauge fields 
and the quarks. First results from simulations 
with an improved fermionic action will however 
be discussed below (for more details see Ref. 0). 

The values of k, L and the (approximate) num- 
ber of configurations used for the calculation of 
the nucleon matrix elements arc the following : 



K 

0.1515 

0.153 

0.155 

0.155 

0.1558 

0.1563 



L 
16 
16 
16 

24 
24 
24 



configs. 
400 
600 
900 
100 
100 
100 



Comparison of the L = 16 and L = 24 lattices at 
k = 0.155 will give us some information on finite 
size effects. 

On the larger lattice we can use lighter quark 
masses without running into problems from se- 
vere finite size effects. This makes the extrapola- 
tion to the chiral limit more reliable, but the time 
needed to invert the fermion matrix increases 
considerably. So it becomes important that we 
choose our algorithm carefully. Therefore we have 
made a comparison between two popular choices, 
minimal residue with overrelaxation (MR) and a 
stabilised variant of the biconjugate gradient al- 
gorithm (BiCGstab). 

Our main findings (for (3 = 6.0, Wilson ac- 
tion quarks) are that if the quarks are not very 
light the MR algorithm, with an overrelaxation 
parameter u = 1.1 was the most efficient al- 
gorithm, saving about 15 % in CPU time com- 
pared with BiCGstab. However as we approach 
k c BiCGstab becomes the preferred algorithm. 



The k where the two inversion times cross over de- 
pends strongly on the lattice size. On a 12 3 x 16 
lattice the crossover is near k ss 0.154 while on 
the 24 3 x 32 lattice BiCGstab does not become 
the better algorithm until k ps 0.1555. For our 
calculations we have therefore used the MR algo- 
rithm for the lower k values, the BiCGstab for 
the highest. 

In order to suppress the unwanted excited 
states as much as possible it is important to 
choose the hadron operators judiciously. We do 
this by applying Jacobi smearing to the stan- 
dard local operators. Both source and sink are 
smeared, since one needs a good projection on 
the ground state on both sides of the inserted op- 
erator in the three-point function. In the case 
of the proton we additionally use the "nonrela- 
tivistic projection" Q. In this way we obtain a 
sufficiently long interval in i, where the two-point 
function is dominated by the ground state. 

Our final choice of the operators whose ma- 
trix elements are calculated is motivated by the 
wish to avoid mixing (as far as possible) as well 
as momenta with more than one nonzero com- 
ponent. Taking the nucleon polarisation (where 
needed) in 2-direction and choosing the momenta 
p = (0, 0, 0), (2n/L, 0, 0) we studied the following 
operators and reduced matrix elements. 



V2,a 
V2,b 
V3 

a 
d 2 
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The two operators for v 2 belong to different rep- 
resentations of H(4). Therefore their comparison 
gives an indication of the size of lattice artifacts. 

Recently, first results have been obtained for 
nucleon matrix elements of operators of the type 



(16) 



which are related to the so-called transversity dis- 
tribution hi pQ|. Purely gluonic operators have 
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Figure 1. R for the pion matrix element V2,b at 
k = 0.153 as a function of r/a. 

also been studied although they fluctuate much 
more than the 2-quark operators 

As explained above, we extract the hadronic 
matrix elements of our operators from the ratios 
R (see (|l3|)). In the case of the 2-quark oper- 
ators we fix i, which allows us to compute the 
quark- line connected part of the three-point func- 
tion for all values of r and arbitrary operators O 
from quark propagators on two sources. The dis- 
connected insertions, which would contribute in 
the flavour-singlet sector only, are however omit- 
ted. We choose t as large as possible in order to 
have enough space for a "plateau" between and 
t where R is independent of r. For the nucleon 
we took t/a = 13 whereas in the case of the pion 
we adopted the symmetrical choice t/a — 16 on 
our 16 3 x 32 lattices. 

Indeed we observe reasonable plateaus. Exam- 
ples for the nucleon are shown in Ref. Q, an 
example for the pion is plotted in Fig.[j} From 
the values of R on the plateaus we can then cal- 
culate the reduced matrix elements taking into 
account kinematical factors, renormalisation con- 
stants etc. 0. 

4. RENORMALISATION 

Using 1-loop lattice perturbation theory we 
have calculated the matrix elements of all op- 




(qu)**2 



Figure 2. Zx Wilson coefficient for the matrix 
element V2, a at K = 0.153. 

erators 0$...^ (n < 4) and 0*$... Mb (n < 2) 
between quark states in the quenched approxi- 
mation |^2| . From these one can immediately 
compute the renormalisation constants and mix- 
ing coefficients for any linear combination that 
one wants to study. Increasing the accuracy of 
the required numerical integrations we obtained 
up to 8 significant digits. In the cases considered 
previously |(|, our results agree with the older 
ones. In the following, we shall use the perturba- 
tive Z's at the scale [i 2 = a~ 2 w 5GeV 2 . 

However, a nonperturbative determination is 
possible ]l3| and should eventually be preferred. 
In the end, the /i dependence of Z must be can- 
celled by the [i dependence of the Wilson co- 
efficient leading to renormalisation prescription 
independent results for the structure functions. 
In order to show to which extent this match- 
ing can be achieved we plot in Figj^ the prod- 
uct of the nonperturbatively calculated Z for V2, a 
with the corresponding renormalisation group im- 
proved nonsinglet Wilson coefficient versus a 2 fi 2 . 
Indeed we obtain a reasonably flat region around 
a 2 /i 2 w 2, though with a rather low value of 
Aqcd ~ lOOMeV. 

5. RESULTS FOR UNPOLARISED NU- 
CLEON STRUCTURE FUNCTIONS 

Plotting our results for i>2, «3, «4 versus 1/k we 
can extrapolate (linearly) to the chiral limit. This 
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Figure 3. (x) for the proton (MOM scheme). The 
circles (boxes) correspond to V2, a (v2,b)- The up- 
per (lower) band of data represents the results for 
the up (down)-quark distribution. 

is shown in Fig.^jfor V2 = (x) in the proton. (Sim- 
ilar pictures for v% and V4 can be found in Ref. §.) 
We see that the values obtained at smaller quark 
masses on the larger lattice (24 3 x 32) are consis- 
tent with those coming from the 16 3 x 32 lattice. 
In particular, at k = 0.155 the numbers from both 
volumes agree within the errors. Furthermore we 
observe approximate consistency among V2, a and 
V2f> indicating the absence of large lattice arti- 
facts (at least in this case). 

At this point it is impossible to resist the temp- 
tation to compare our results with experimen- 
tal numbers. However, due to our use of the 
quenched approximation we can only compare 
with valence quark distributions. But even these 
will be influenced by the presence (or absence) 
of the sea. In fact, the phenomenological values 
indicated by the asterisks in Fig.|| differ substan- 
tially from the extrapolated lattice data. One 
might however expect that the flavour-nonsinglet 
combination {x}^ — is less sensitive to sea 

quark effects. Indeed our result of 0.23(3) (MS 



scheme) obtained by averaging v 2 , a and v 2 ,b com- 
pares more favourably with the phenomenological 
number 0.18. 

On the other hand, the large discrepancy be- 
tween the quenched and the phenomenological 
values of (x)^ should not be too surprising. 
In the real proton the contributions of the va- 
lence quarks, the sea quarks, and the gluons 
must add up to 1. In the quenched approxima- 
tion, the missing sea contribution, which is about 
0.18, must somehow be compensated by the va- 
lence quarks and the gluons. Assuming that the 
gluon distribution is not too strongly affected by 
quenching one would expect the quenched (x)^ 
to be larger than the phenomenological result, 
which is exactly what we find. 

The values obtained for (x 2 ) and (x 3 ) (MS 
scheme) differ less drastically from their phe- 
nomenological counterparts: 



lattice 



experiment 



(x)a 

(x) b 

<* 3 > 



0.369(26) 
0.440(22) 
0.108(16) 
0.020(10) 



d 

0.169(13) 
0.187(10) 
0.036(8) 
-0.001(6) 



u 

0.284 
0.284 
0.083 
0.032 



d 

0.102 
0.102 
0.025 
0.008 



6. RESULTS FOR POLARISED NU- 
CLEON STRUCTURE FUNCTIONS 

As in the case of the unpolarised structure func- 
tions, the chiral extrapolation of the 16 3 x 32 data 
is confirmed by the results obtained on the 
24 3 x 32 lattice at smaller quark masses. This is 
shown in Fig.[| for the case of Aq. In the chiral 
limit wc find 



Am = 0.84(5) , Ad = -0.24(1) . 



(17) 



These numbers should be compared with the 
fraction of the proton spin carried by the valence 
quarks (see, e.g., Ref. [H): 

Au v = 0.92 , Ad v = -0.34 . (18) 

Again one might hope that the difference of u and 
d contributions is less sensitive to quenching. One 
finds Am — Ad = 1.08(6) to be compared with the 
axial vector coupling constant gA = 1.26. 

In the higher moments of g\ sea quark effects 
are also expected to be suppressed so that the 
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Figure 4. Am (upper values) and Ad (lower val- 
ues) for the proton. 



Figure 5. 



) for the pion versus 



{m^/mpY for n = 2,3,4 (from top to bottom). 



quenched results should be reasonably close to the 
experimental numbers. For the proton we find 

/ dxx 2 g 1 (x 1 Q 2 ) = 0.0150(32) , (19) 

JQ 

where Q 2 = [i 2 « 5GeV 2 . The E143 collabora- 
tion obtains in a recent analysis 0.0121(10) 

7. PION STRUCTURE FUNCTIONS 

In FigjH] we summarise our results for the va- 
lence quark distribution in the pion obtained on 
a 16 3 x 32 lattice |T(|. We used the same oper- 
ators as in the nucleon case. Like for (x)^ in 
the proton, v 2y b is consistently larger than Vi a . 
The three filled squares on the right represent 
the heavy quark limit, whereas those on the left 
are calculated from a phenomcnological valence 
quark distribution Jl7[ . 

8. FIRST RESULTS FROM AN IM- 
PROVED ACTION 

Wilson's fermion action suffers from 0(a) lat- 
tice artifacts. In order to suppress them at least in 



on-shell quantities, Sheikholeslami and Wohlert 
|l8|| added the so-called clover term 

-cs W ay^ j q(x)cx IJ , v F IJ , v (x)q{x) (20) 

to the action, where is the clover-leaf lattice 
version of the field strength. For the coefficient 
csw one obtains in lowest-order perturbation the- 
ory c sw = 1- A nonperturbative determination 
is however possible as the ALPHA collaboration 
has shown fl9| . For (3 — 6.0 they find the optimal 
value csw = 1-769 @, which should reduce the 
discretisation errors from O(a) to 0(a 2 ). Since 
a "canonical" value of csw has not yet emerged, 
we have extended our perturbative calculation of 
renormalisation constants for the local 2-quark 
operators to the Sheikholeslami- Wohlert action 
with arbitrary csw ||l| • 

In order to improve matrix elements such as 
those considered here one needs in addition im- 
proved operators. We have obtained first re- 
sults with perturbatively improved operators on 
a 16 3 x 32 lattice at (3 = 6.0 and c sw = 1-769 ||. 
The value of Au— Ad in the chiral limit changes to 
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1.22(14) for csw — 1-769 leading to better agree- 
ment with the experimental result 1.26. 

9. CONCLUSIONS 

This talk has described recent attempts to cal- 
culate hadron structure functions on the lattice. 
In the quenched approximation, at least the lower 
moments can be calculated with reasonable sta- 
tistical accuracy. 

Concerning the systematic uncertainties we can 
make the following statements: 

• The comparison of results from 16 3 x 32 and 
24 3 x 32 lattices does not reveal large finite 
size effects. 

• The extrapolation to the chiral limit seems 
to be smooth. In the range of quark masses 
that we studied we did not observe any un- 
expected mass dependence with the possi- 
ble exception of c?2 • 

• As a check for cut-off effects we calculated 
(x) from two different lattice operators. 
The outcome indicates that these effects are 
not too large. First results obtained with a 
nonperturbatively improved fermion action 
look promising. 

Still there are further sources of systematic er- 
rors due to the renormalisation constants, the 
contributions of purely gluonic operators and 
fermion-line disconnected parts of 2-quark oper- 
ators, etc. A major problem is, of course, the 
quenched approximation. However, the overall 
agreement with the real world is already rather 
satisfactory, at least as far as quantities are con- 
cerned which are expected to be less sensitive to 
quenching. 
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